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Abstract 

In this work we incorporate some real-system effects into the theory of ori- 
entational phase transitions under shear flow (M. E. Cates and S. T. Milner, 
Phys. Rev. Lett. 62 1856 (1989) and G. H. Fredrickson, J. Rheol. 38, 1045 
(1994)). In particular, we study the influence of the shear-cell boundaries on 
the orientation of the lamellar phase. We predict that at low shear rates the 
parallel orientation appears to be stable. We show that there is a critical value 
of the shear rate at which the parallel orientation loses its stability and the 
perpendicular one appears immediately below the spinodal. We associate this 
transition with a crossover from the fluctuation to the mean-field behaviour. 
At lower temperatures the stability of the parallel orientation is restored. We 
find that the region of stability of the perpendicular orientation rapidly de- 
creases as shear rate increases. This behaviour might be misinterpreted as an 
additional perpendicular to parallel transition recently discussed in literature. 
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I. INTRODUCTION 



When subjected to a shear flow, AB block copolymer melts exhibit an orientational phase 
behaviour which is absent in equilibrium. A system under shear shows not only transitions 
between different morphologies (typically lamellar, hexagonal, cubic and gyroid but 
also transitions between different orientations of these morphologies with respect to the 
shear geometry Experimental literature extensively discusses this effect for lamellar ^SM 
and hexagonal phases 

The theoretical description of the lamellar reorientation was developed in fftffl. The same 
method was applied in [[| to study the hexagonal pattern. In these theories orientational 
transitions appear as a result of interaction of shear flow with critical fluctuations in melt. 
There are two distinct regimes: a slow flow only slightly perturbs the fluctuation spectrum 
while a fast flow significantly dumps fluctuations, restoring the mean-field behaviour in the 
limit of infinite shear rate D — > oo. Correspondingly, the parallel lamellae (their normal 
is parallel to the shear gradient direction) are found to be stable in the small shear rate 
regime, while the perpendicular lamellae (their normal is perpendicular to both the gradient 
and flow directions) are stable at high shear rates. Fredrickson has shown that if one takes 
into account the difference in viscosities of the pure melt components, the perpendicular 
phase loses its stability at low enough temperatures and the parallel orientation is restored. 
Schematically this behaviour is summarized in Fig.[l| 

However, there is an experimental evidence that this picture is not complete. At very 
high shear rates the parallel orientation was found to be the only stable one 0,0. This 
cannot be explained in the discussed framework of ||, since it predicts the stability region 
for the perpendicular phase to increase as D — > oo. 

In this work we propose an explanation of the additional transition (C-transition in 
Fig.|I|). We argue that the missing element of the theory is the interaction of the block 
copolymer melt with the walls of the shear cell. We consider a block copolymer film con- 
fined in-between two walls in the gradient direction and subjected to a steady shear flow. 
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Usually the distance between the interfaces in the other two directions is much larger and 
we ignore their influence. This model will predict the parallel orientation to be stable in 
the D — > oo limit since the influence of shear and fluctuations vanishes in this limit. The 
only symmetry-breaking factor is then the wall-copolymer interaction which stabilizes the 
parallel orientation [|T^,|T^]. The complex behaviour at lower shear rates will arise from the 



interplay of three factors: shear flow, fluctuations and wall-melt interactions. 

We admit that the influence of the surface interactions is possibly small. However, 
Balsara et. al reported jl4| that in the absence of shear the walls of their shear cell induced 



the parallel alignment through the whole 0.5-mm-sample, although the lamellar spacing is 
somewhat 4 orders of magnitude smaller. Under shear Laurer et. al |15| observed that 
independently of the bulk orientation there is always a near-surface layer of the parallel 
lamellae which penetrates up to 2 /zm into the bulk. Thus, even a weak symmetry-breaking 
field can be crucial in the absence of other factors. 

We also want to mention that the equilibrium theory of block-copolymer melt ordering 
near surfaces is well-developed fll2] , |l3l , |rB| -|2"2]] . Some questions about dynamics of such an 
ordering were addressed in P3"| , |2T)|| . However, until now this theory was never applied to 
non-equilibrium systems. 

Our paper is organized as follows. In Section II we derive the equations governing the 
dynamics of the melt and construct a non-equilibrium potential whose minimal value will 
determine the stable orientation. In the first part of Section III we estimate the shear rate of 
the A-transition while the other two transitions (B and C) are analyzed in the second part. 
In conclusion we discuss in detail properties of the obtained phase diagram. In Appendix 
we provide an example clarifying the role of thermal fluctuations. 

II. DYNAMIC EQUATIONS 

Let us consider a block copolymer melt confined in-between two surfaces in the y- 
direction. It is also subjected to a steady shear flow v = Dye x (see Fig.^). We ignore 
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any alteration of this velocity profile and assume that it is kept through the whole system. 
We choose the local deviation of composition from its average to be an order parameter 0(r) 
and define its Fourier transform as 



where 



0(k) 



dr e 



— ikr , 



and 



e ikr 0(k) 



(2.1) 



d kg- dkg 
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dr 
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oo J-L/2 J -oo 



(2.2) 



It is convenient to work in dimensionless units and we rescale lengths and wave-vectors: 
r — > 6 _1 r and k — > 6 k, b being the size of a monomer. 

Following [^4|,0,|[ we assume that the dynamics under shear flow is governed by the 
Fokker-Planck equation: 

5 5H 



dP[ M _ f 5 



Of 



5<j){r) 



5<p{r) 



r) J ' ^ dx 



P[r,t] 



(2.3) 



where P[4>,t] is the probability to realize the order parameter profile 0(r) at time t, and fi 
is an Onsager coefficient. In Appendix we provide arguments for using the Fokker-Planck 
equation instead of any other deterministic equation. 

In eq. (|2.3| ) the Hamiltonian TC consist of two contributions: the bulk Hamiltonian derived 
by Leibler [§ 

XH L [<t>] = \ /r 2 (g)0(q)0(-q) + ^ I I I r 3 (qi,q 2 ,q 3 )0(qi)0(q 2 )0(q3) 

^ Jq J qi Jq 2 J qz 

T\ [ [ [ [ r 4(qi ; q 2 ,q3,q4)0(qi)0(q 2 )0(q3)0(q4) (2.4) 



+ 



and the surface energy [|12| , [13 1 
NH< = [dr 



fl^(r) + |0(r) 2 



(2.5) 



where iV is a number of monomers in a molecule, i?i ~ (x^)cop-surf is the strength of 
the interaction between the surface and the copolymer melt and describes the additional 
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interaction in the melt induced by the presence of the surface (it changes the local tem- 
perature in the vicinity of the surface). Our goal is to construct a real-space version of the 
Leibler Hamiltonian Hl- In p5| , |26|JT^| it was shown how to deal with the second-order vertex 
function. Separating small- and large-wave-vector asymptotic behaviour, one can show that 

A 



T 2 (q)^^ + Bq 2 -x (2.6) 



where 



.4 



2RlP{l-fY 



B = v^T) (2 ' 7) 



/ 3 (1-/) 3 



1/2 



with / = Na/N being the volume fraction of the A-component. In |L2| the third- and 
fourth-order vertex functions were assumed to be constant. However, as it was noticed in 
|8|j9|j2^], it is crucial to keep the angle-dependence of the fourth-order vertex function in 
order to discriminate between the parallel and the perpendicular orientations. There, the 
following approximation was made: 

iMoi, q2, q3) = 5(qi + q2 + q3)r 3 

r 4 (k,q,-k,-q) = A(l-/?(k-q) 2 ) , (2.8) 



where k = k/fc. In eq.( |2.8| ) all the wave- vectors are assumed to have the same length 
go = y A/B, which corresponds to the first unstable mode on the spinodal 0. The assump- 
tion /3 < 1 was shown to be correct for almost every architecture of AB block-copolymer 
molecules [^] (for example, for diblocks (3 < 0.1). For an arbitrary star of 4 q's one can 



write to the lowest order in angles [28,29 



r^qi, q2, q3, q4) = 5 (qi + q2 + q3 + q4) x 
Ao + Ai ((qi • q 2 )(qs ■ q4) + (qi ■ qs)(q2 ■ q4) + (qi ■ q4)(q2 ■ qa)) (2.9) 



+A 2 ((qi • q 2 ) 2 + (qi • q 3 ) 2 + (qi ■ q 4 ) 2 + (q 2 • q 3 ) 2 + (q 2 • q 4 ) 2 + (q 3 • q 4 ) 



Ai A 2 
Y ' Y <<c 1 



Comparison with eq. (|2.8| ) gives 
A = A + Ai + 2A 2 



2Ai + 4A 2 



A + Ai + 2A 2 

Thus, the required real-space representation of the Hamiltonian Tt can be written as 

NH[<j>] = Jdr [| (V0(r)) 2 - I^(r) a + | |«(r - r')0(r)0(r') + ^(r) 3 + ^0(r) 
3Ai (v0(r)-V0(r)) 2 6A2 
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where 



S(r-r') 



3 iq(r-r') 



h(r)(j)(r) 



+ NH S (2.11) 



(2.12) 



Here we have added an auxiliary field h which will help us to construct a thermodynamical 
potential governing the dynamics under shear. Afterwards it will be set to zero. 

The Fokker-Planck equation (|2.3|) together with eq ns . ( pTTT] . [2 . T| . [2~5|) form a phenomeno- 
logical set of equations describing the dynamics of block-copolymer melt under shear flow 
in the presence of surfaces. We do not solve these equations directly, but following M we 



use the method of Zwanzig |3(J to derive a system of coupled equations for the first two 
cumulants of P[4>, t] 



c r 



(0(r)> 

5(r 1 ,r 2 ) = (0(r 1 )0(r 2 ))-(0(r 1 ))(0(r 2 )) 



(2.13) 



where c is the average order parameter profile, and the structure factor S is a measure of 
the fluctuation' strength. We introduce a generating functional 



G[£,t] = log y^exp 



dr<f>(r)£(r) 



P[(f>,t] , 



(2-14) 



use eg. (|2~3|) to derive an equation of motion for and then expand this equation in 

terms of £. The two lowest-order equations read: 



1 (M r ; --y^- + BAc ( r ) + Xc(r) - A /dr' 0(r - r')c(r') - — 
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c(r) V Q V^c(r) + c(r)5 + c(r) 2 V a V^c(r) + 5(0)A 2 c(r 



(2.15) 



V a V^c(r) 5 a/3 + h{r) + [H x - a lC (r 
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1 <95(r — rn) . D 5(r — ri) _ . . _,. . 
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5(r-n) V a V^c(r) (2.16) 



+ ( V a 5(r - n) J (Vc(r)j + 5^V Q 5(r - r 1;j ^ , 

+5(r - n)5 + 2c(r) (v Q V^c(r)) (v„V^(r - n)) + 5(0) A 2 5(r - n 
+2V ct V /3 (c(r)5(r - ri)V a V/9c(r)) + V a V^ (c(r) 2 V a V^5(r - n 
-a 1 S(r-r 1 ) S (y + - J + 5 (y - ^ 



where 



£tf = V^V' 5(r - r') 



S = V Q V^V^V^5(r - r') 



(2.17) 



Here we have neglected all higher cumulants and made use of a natural assumption 
S(r 1 ,r 2 ) = S(ri - r 2 ). 

Apart from the surface terms, eqns.( |2TT5]j2.16| ) are the real-space analog of the eqns.(2.25- 
26) from ||. Here the terms proportional to 5(0) play the role of the fluctuation integral 
<x(k) from §: 



a(k) = ^Js(q) l-/3(k-q) 



(2.18) 



To keep our model as simple as possible we leave only the linear term in the surface 
energy ( |2.5| ) and put a x = 0. Then we set 

c(r) = 2a cos(g n • r + <p) (2-19) 
h(r) = 2h cos(g n ■ r + <p) 

where a is yet to be determined amplitude, n is a unit vector perpendicular to the surface 
of the lamellae and (p is a phase shift which will be chosen to minimize the surface energy. 
The auxiliary field h simply follows the behaviour of c. Fredrickson has shown |12j that in 
equilibrium the presence of the surfaces causes spatial variations of the amplitude a which 
decay exponentially away from the surface. Since we are only interested in the orientation of 
the lamellar profile ( |2.19| ), we ignore the spacial dependence of a and set it constant. With 
these simplifications the equation for the Fourier transform S(k) of S(ri — ^) from eq. fl2.16Q 
reads 

2/i dt 2/i dk y v ' 



where 



S \k) = r - k-e -k + Bk 2 + ^ - Xs 

k e k = 2(( X JV), - XN) + Xa 2 (l - (3(n ■ k) 2 ) + <r(k) (2.21) 



Here we have introduced the same notation as in 0|§. In eq.( |2.21| ) S'o(k) is the equilib 



rium structure factor and r — k- e -k denotes the renormalized temperature. Within the 
fluctuation theory the spinodal temperature determined from the condition 



r — k- e k 



a=0 



(2.22) 



differs from the mean-field value {x^j ■ In the case (3 = such a fluctuation correction was 
discussed in |JT| . The presence of shear breaks the rotational symmetry and the spinodal 
temperature becomes orientation- dependent. This gives rise to the — k- e -k term, with 
e ij ~ P ( see eq-(|2-(j|))- Here the role of the angle- dependency in T 4 is especially transparent: 



if /3 = 0, we would not be able to discriminate between different orientations. 



The method of characteristics [32] gives a formal solution for the eq.( [2.20|) : 
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(2.23) 



The steady-state regime is approached as t — > oo. The integration in Q2.23 ) can be performed 
in the limiting cases D —>■ and D — > oo and will be discussed in the next section. 

Now we derive an equation for the amplitude a. We substitute the lamellar profile Q2.19|) 
into eq.( 2.15 ) and perform an averaging over the lamellar period 



n x n y n z q% / v 
( --' ) = ^2^ / <IJ ' 



2tt 



2tt 
<3Q n z 



dzcos (gon • r + <£>) 



(2.24) 



Discarding the transverse orientations with n x ^ |p4|,[7|,|8| , we obtain 



— — = h — (r — n- e -n)a + -A(l — (3)a 3 + rjcos(ip)8, 



n|,l 



(2.25) 



where r\ = ^Hi, and 5 n 2 ;1 is the Kronecker delta-symbol which is non-zero only for the 
parallel | = 1) orientation. Following || we notice that the equation ( |2.25| ) has a 
gradient form (with h = 0): 



1 da 1 9$ 
// dt 2 da 

Since the potential <J> can only decrease with time: 

<9$ 9$ 9a « / (9$ 



(2.26) 



<9t (9a 9t 2 \ da 



< 



(2.27) 



the steady-state of the system will be determined by the minimum of $. Now we use the 
auxiliary field h to construct $. In steady-state da/dt = 0, and $ is obtained by integrating 

Id* 



h 



2 da 



(2.28) 



Using h from eq.( 2.25|) , we obtain 



$ = $ - 2r]aS n 2 x 



(2.29) 



where 



1 /■"<-» 
$ = --A(l -/3)a 4 + 2 / da'(r - n- e -n)a' (2.30) 

In eg. ( P-29| ) we have already minimized with respect to the phase shift <p, assuming that 
a > (the other terms depend only on even powers of a and are not influenced by this 
choice) . 

The non-trivial dependency of r — n- e n on a comes from the term proportional to 
cr(k) in ( [2.21| ) and the potential $ appears to be dependent on the fluctuation integral via 



eq.( |2.30| ). Now we are ready to discuss the stable orientations in different regimes. 

III. PHASE TRANSITIONS 
Crossover from small- to high-shear rate behaviour 

In this subsection we analyze the transition from the parallel to perpendicular orienta- 
tion caused by increase of shear rate (the A-transition in Fig.|l|). We start with noticing 
that at low shear rates the parallel orientation is the only stable one. Indeed, as it was 
shown by Fredrickson ||, <£>o is minimal for n 2 = 1 in the limit D — > 0. The surface 
term in eq. ( |2.29| ) also favours the parallel orientation. Thus, our theory does not modify 
Fredrickson's prediction for small shear rates. 

At high shear rates D — > oo, the integration in eq. ([2.23|) can be performed |24j], yielding 



where 



2 \ 2 /3 

s ~< k >-°i7^y (3 - i) 



Co= (075 and Q = 1r (3 - 2) 



For the intermediate shear rates 5(k) can be interpolated between Sq and S 1 , 



S(k) 



- <-» - A 
r-k- e -k + Bk 2 + — - Xs + 

k z 



l_ ( yf0LD\k x ky\ \ 

co V ml ) 



2/3" 1 



(3.3) 



One should realize that the previous equation is an analytic continuation of the D — > oo 
behaviour to D < oo values. As a result, a small-/} behaviour of eq. (|3.3|) does not correspond 
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to the D — > behaviour of eq. (|2.23|) . On contrary, it describes the D ~ 0(1) region. Since 
we expect the A-transition to lay in-between the D ~ 0(1) and D — > oo regions, we need to 
calculate the fluctuation integral er(k) in-between these regions. This can be done in several 
steps. First, we use 5(k) from eq.( |3.3| ) to perform the radial part of the integral in eq.( [2.18| ). 
This gives 



V ; 16vr 3 



1/3 



1-/3 (k-q 

It £ 1 1/3 



7T 



— arct.au f ( — ^= I 

2 I IA;^|V3 \D^J 



1/3" 



(3.4) 



As a next step we expand the integrand for D <C 1 and D ^> 1 and sum these expressions 
keeping only the few first terms. Integration over the orientations of the unit vector q 
(J dtt = d9 sin 9 f^dtfi) then gives 

A /n f / ?\ 1 



+4vr 2 2 1/3 \/3Z 







1-^(2^ + 2^ 



3A; 



3Z 3 L 



/ 1 -/5(/ 2 (^ + ^)+/ 3 A; z 2 
h-/3(h(k 2 x + k 2 v ) + I 6 k 2 z 



(3.5) 



where 



V47TCQ 



h 



, y^r(|) : 
^ r(f) 
_r(|) 3 



.1/3 

A D\ 

2 ^r(|)r(f) 

r(f) 
v^r(i)r(|) 



D* = \fiyfa 



„ r(l) 2 r(l) 

r(f) 
r(l)r(i) 2 

' r(¥) 



r(¥) 

This procedure is ill-defined from mathematical point of view. However, any direct analysis 
of eq.( |2.23| ) is impossible and we use eq.( |3.5| ) for moderate shear rates. 

Since eq.( |3.5| ) does not depend on the renormalized temperature, the integration in 
eq.( |2.30| ) is trivial and gives 



$ = [r + a(n)] a 2 + -Aa 4 (l - (3) - 2 V a5 n%il 



(3.6) 



where 
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Minimization with respect to a gives to the first order in r\ 

_ [r + a(n)] 2 / o * + 

$ " " A(l-/3) " 2T] \j " 2 W^T) 5 < l (3 ' 7) 

The order-disorder transition (ODT) occurs when <3> becomes negative. The corresponding 

transition temperature is 

r s (n) = -<r(n) (3.8) 



which coincides with eq.( |2.22| ). The orientation with the lowest a will appear immediately 



below the ODT temperature. The crossover (the A-transition) is then located at such a 
value of D that cry — cx± changes its sign. From eq.( ^.5[ ) this point is given by 



cry - a± 



3^(/ 3 - h) + ^2^V3Z + Z\h - h) = (3.9) 



and is found to be 



D cr « 2.64/za (3.10) 

From eq.( |3.9| ) it also follows that rj > for D < D cr , which fits the small-shear behaviour 
discussed in the beginning of this subsection. When D > D cr , the perpendicular orientation 
first appear below the spinodal. This crossover is depicted in Fig.|3|. 

As it was noticed before the mean-field behaviour is restored in the limit D — > oo. 
On the other hand, the small-D region is dominated by fluctuations. Thus, D cr can be 
interpreted as a position of a crossover from the fluctuation to mean-field behaviour. The 
scaling properties of D cr follow from eq. (|3.10|) and are determined by the Onsager coefficient 
/i. Using the results of |33|-|36| (/i = qoX(qo)/N with A from [ |34l , j35[| ) we obtain 



D cr ~N- 3 (3.11) 

which shows that the fluctuation region disappears in the limit N —>■ oo. In equilibrium the 
same conclusion was drawn in [p,|3T|. 

Finally, we emphasize that the results of this subsection are independent of the surface 
interaction. The same results can be obtained within the Fredrickson theory |§ (rj = 0). 
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B- and C-transitions 



In the previous subsection we have discussed the order-disorder transitions. Now we 
consider lower temperatures and look for transitions between different orientations in the 
high-shear limit. The corresponding free energies are given by eq. ([3/ 

, 2 



A(l-/3) 



A(l-/3) 



(3.12) 



r + a ± 



A(l-/3) 

where cx(n) is given by its high-shear limit of eq.(|3 

,A) 2/3 / D*\ 1/3 2V3V3 



CI/ 



a n 



v ~ 1-^(2^ + 3^) 



(3.13) 



To the leading order in D*/D, the transition from the perpendicular to parallel orienta- 
tion occurs at temperatures which are the roots of the equation $y = $^ 

^ 4 2r/ 2 A(l-/5) 



Ti = -an 



T"2 



(<7|| - CX ± 



if r) 2 \(l-(3) > (a,, -a ± ) 3 (3.14) 



where 



(Til — <7j_ = (3 



2 1 /3 V ^ ( a A) 2/3 /AA 1/3 



56 



(3.15) 



There ri corresponds to the I — 1 1 transition, while r 2 - to the reverse one. Now we summarize 

our results in a phase diagram. 



IV. DISCUSSION OF THE PHASE DIAGRAM AND CONCLUSION 

In this work we incorporated some real-system properties into the previously developed 
theory of the orientational phase transitions under shear flow 0^0- In particular, we consid- 
ered the influence of the shear-cell boundaries in the gradient direction on the orientation of 
the lamellar phase. In equilibrium the lamellae are known to orient parallel with respect to 
the boundaries |12 , |13| . Under shear the tendency to orient parallel to the surfaces competes 
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with the orientation favoured by the flow which appears as a result of the coupling between 
the flow velocity field and the order-parameter fluctuations The interplay between 

these two factors produces the non-trivial phase diagram shown in Fig.|j. 

At low shear rates the parallel orientation is preferred by both the shear and surface 
terms in eq.( |2.29| ). Therefore it is the only stable orientation in that part of the phase 
diagram. When shear rate reaches the value D cr given by eq. (|3.10|) , the perpendicular 
orientation becomes stable immediately below the ODT temperature. We associate this 
change in orientation with the crossover from the fluctuation-dominated behaviour to the 
mean-field one. Indeed, at very small shear rates the equilibrium fluctuation spectrum is 
only slightly modified by the flow, while at high shear rates the flow strongly suppresses 
fluctuations and restores the mean-field behaviour. Therefore, there is a crossover point and 
the corresponding change of orientation. 

At high shear rates and away from the spinodal, the surface influence starts to play an 
important role. In the narrow region between D cr and Di, estimated from the condition in 

eq.(EH) 

Di = 2 (3<^^ Ac* 



56 3 1 -0 r] 2 B 9 / 2 

the influence of shear is still very strong and is capable of stabilizing the perpendicular 
orientation at all temperatures. The size of this region is very small due to the scaling 
D 1 ~ iV- 13 / 2 . When D > D 1 there appears a region where the parallel orientation is stable. 
It takes over the perpendicular one at r = T\ and looses its stability again at r = r 2 given 
by eq.( |3.14 ). This region grows as the shear rate increases, and in the limit D = oo the 
parallel orientation occupies the whole range of temperatures (0, — oo). This coincides with 
the predictions of the equilibrium mean-field theory [[l2i[T3| . We therefore argue that there 
is no sharp C-transition as shown in Fig.|I[ Since the region between the spinodal and the 
parallel phase shrinks with an increase of the shear rate, there always be some value D 2 such 
that for D > D 2 the size of this region will be smaller than the resolution of the experimental 
device. This value D 2 can be misinterpreted as a position of an additional transition. 
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An important feature of our theory is that it is able to reproduce the B-transition without 
additional assumptions. In the previous theory ||, Fredrickson had to take into account 
the difference in viscosities of the pure components in order to reproduce the B-transition. 
Namely, he put rj [4>] = ?7o+^i0 which can be considered as a Taylor expansion of the viscosity 
77 [0]. As a result, in high-shear limit the size of the stability region for the perpendicular 
phase is of order of {t]o/vi) 2 an d grows as D — > 00. While depicting the main physics, 
this approach has internal problems since the derivative 7/1 is not a well-defined object and 
therefore the whole theory depends on a phenomenological parameter which is difficult to 
estimate. Moreover, Fredrickson's theory does not predict the C-transition. Our theory is 
free from these problems. It, however, predicts the ||— >_L transition at very low temperatures 
in high-shear limit which was not observed. There could be several explanations of this 
prediction. First, this transition occurs at very low temperatures (r 2 ) where the weak- 
segregation theory does not work. Second, this transition might be an artifact of the 0(r]) 
expansion. Finally, this transition can be removed if we use Fredrickson's argument about 
the viscosity dependence on the order parameter. It stabilizes the parallel orientation at low 
temperatures. 

It is possible that the absolute value of the surface interaction is small. However, since 
it acts as a symmetry-breaking factor, its influence is very important |T^Jl~5[| . This statement 



can be checked experimentally. We have shown that the positions of the transition lines 
in the high-shear limit are dependent on the strength of the surface-copolymer interaction 
7] ~ (xN)cop-surf- Therefore, the phase diagram of a particular copolymer system depends 
on the material of the shear-cell walls. A systematic study of this dependency will provide 
arguments for or against our theory. 

In this work we have considered the influence of the walls in the gradient direction. We 
also want to comment on the role of the boundaries in the other shear directions (flow and 
vorticity). Formally, these walls will also induce alignment parallel to themselves. However, 
the flow profile near those walls is no longer a simple triangular one and we expect this 
disordered flow to destroy their orientational tendency. Moreover, the distance between 
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those surfaces is normally much larger than between the walls in the gradient direction and 
their influence is thus weaker. Therefore we neglected them in our work. 

Finally, we want to discuss briefly possible modifications and extensions of the developed 
theory. A very interesting problem is to calculate the alterations of the density profile in a 
confined system under shear. This can be achieved by restoring the position dependency in 
the amplitude a and deriving the corresponding amplitude equation from eq.( [2.15|) . In the 
absence of shear this problem was solved in [12| . Another possibility is to use our formalism 



for other external fields rather than interactions with surfaces. A good example is an electric 
field which is coupled to the square of the order parameter [|37],[^| . With some modifications 
eq. (|2.15| , |2.16|) can be a starting point for the corresponding theory. Importance of such a 



theory for a system in electric field and simultaneously under shear was outlined in [38 
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APPENDIX: INFLUENCE OF THERMAL FLUCTUATIONS 

Here we want to emphasize the role of fluctuations in our theory. We are going to show 
that if the Fokker-Planck equation ( |2.3| ) is replaced by a deterministic one, the theory will 
not be able to discriminate between different orientations. What follows should not be 
considered as a proof but, more likely, as an illustration that has general features. 

Let us consider a Langevin equation equivalent to the Fokker-Planck equation fl2.3Q . If 
we now remove the noise term, it reads 

_E + V .V^ = ,V*— (A.1) 

where an Onsager mobility fx differs from /i in eq. Q2.3p and H,l is the Leibler Hamiltonian 
(eq. Q2.llj) ) without 7i s and h = 0. A similar equation was considered in | 39| . There the 
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authors used a Hamiltonian with Ai = A2 = and showed that in steady-state the theory 
predicts both orientations to be equally stable at all shear rates. This is not surprising since 
their theory does not contain fluctuations and the angular dependence of the fourth-order 
vertex function r 4 - the two ingredients that were argued to be crucial in explaining the 
reorientation phenomena PJ5j,P^]. 

In order to separate these two effects we keep the angular dependence in T 4 (Ai, A2 7^ 0), 
but use the deterministic equation (|A.1|) . We follow the approach of |39[ and derive an 
amplitude equation from eq. (|A.l| ) assuming a single plane- wave density profile. The solution 
of this equation is 

drf{r)p-\r) ^ 
'0 J 



A(t) 



P[l) ' 3P(f) fdrf{r)p-\r) 



A(0) 



(A.2) 



and 



/(*) 



4Q 



1 z Ai + 2A2 



An 



P(t) = exp 



a(t) = (l + e)Q 2 -Q 4 



Q 2 = ql + (q y ~ Dtq x f + q\ 



where A is the amplitude, e = (x — Xs)/Xs, and the amplitude and the units of length and 
time are scaled with a/x s /A , ^/B/xs and B/(jlxl), respectively. We see that the equation 
becomes symmetric with respect to the interchange q y q z . Thus, even in the presence 
of the angular dependence in T 4 it is impossible to distinguish between the parallel and 
perpendicular orientations starting from eq.( |A.l| ). 

We believe that a theory without fluctuations of the order parameter is not capable of 
describing the reorient at ional transitions. There are phenomena where fluctuations only 
modify a deterministic behaviour [fl~0|-|4~^l . However, the orientational behaviour under shear 
flow does not belong to this class of phenomena. It can only occur in the presence of thermal 
fluctuations. 
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APPENDIX: FIGURE CAPTIONS 



Fig.l. Schematic phase diagram "temperature vs. shear rate" as a compilation of the 
theoretical predictions by Fredrickson and the experimentally observed C-transition. The 
order-disorder transition line behaves as r ~ D 2 for small shear rates D and r ~ D^ 1 ^ for 
D — > oo. The B-transition line levels off in Fredrickson's theory. 

Fig. 2. Orientations of the lamellar phase in a simple shear flow. The axes of the coordi- 
nate system correspond to the shear geometry: x is the flow direction (v), y is the gradient 
direction (Vv x ) and z is the vorticity direction (V x v). In the parallel orientation the nor- 
mal to the lamellar layers is oriented parallel to the gradient direction, in the perpendicular 
- to the vorticity direction. The walls in y direction interact with the melt and prefer one 
of the components. 

Fig. 3. Schematic behaviour of the spinodal temperatures for the parallel and perpendic- 
ular orientations in the vicinity of the crossover point. 

Fig. 4. Phase diagram for the lamellar phase under steady simple shear flow as predicted 
in this work. 



21 




-X u 



fluctuations mean-field 




